February 1969 


Brief 69-10038 




NASA Tech Briefs are issued to summarize specific innovations derived from the U.S. space program, 
to encourage their commercial application. Copies are available to the public at 15 cents each 
from the Clearinghouse for Federal Scientific and Technical Information, Springfield, Virginia 22151. 



Monte Carlo Direct View Factor and Generalized Radiative 
Heat Transfer Programs 


The problem : 

To develop a method to find the view factor from 
a surface segment to another surface segment and to 
find Radiative-Transfer Coefficients between the sur- 
face segments. 

The solution: 

A package of computer programs was written: (1) 
to find, using the Monte Carlo technique, the “black 
body” view factor from one surface segment to another 
surface segment defined on a primary surface; and (2) 
to define the primary surfaces (collection of simple 
geometric objects) in three-dimensional space, and 
find the radiative-transfer coefficients from a surface 
segment defined on a primary surface to all other sur- 
face segments defined in a real enclosure. 

How it's done: 

The Monte Carlo Direct View Factor Program uses 
the Monte Carlo technique to find the “black body” 
view factor from a surface segment defined on a pri- 
mary surface to another surface segment defined on a 
primary surface. The collection 'of geometric objects 
(primary surfaces) includes a sphere, cylinder, cone, 
disc, and parallelogram. 

The function of the program is numerical evalua- 
tion of the integral: 
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where x ( , Xj are veciors representing points on surfaces 
I, J of area A,. ( . B (“j. Xj) is a blocking factor which 
is unity if X| can ’ % see’* Xj, and zero otherwise. H(xj. 
Xj) = n, • (Xj-x,) nj -(x, -Xj) / tt| x , — x J . where n, and 
fij denote unit surface normal vectors. 


The segment (or node) object numbers are obtained 
from the view factor command. The nodes are located 
in the definition table. From the nodes come the pri- 
mary surface object numbers and K values. The pri- 
mary surfaces are located in the definition table and 
copied into temporary storage. An error message is 
printed if nodes or primary surfaces are missing. The 
areas of the nodes are calculated and the canonical 
forms are generated for the subroutines: 

F Reset the interation sums and parameters 

2. N = N + l 

3. Generate random points and Xj. 

Two random numbers are obtained from a uniform 
distribution. These numbers are modified so that the 
points generated on the different surface window will 
cover the surfaces evenly. The surface points and nor- 
mals are obtained using the subroutines. 

The Generalized Radiative Heat Transfer Program 
finds radiative-transfer coefficients from a surface 
segment defined on a primary surface to all other 
surface segments defined in a real enclosure. The pro- 
gram permits radiation heat transfer analysis for en- 
closures containing surfaces which have, in general, 
both diffuse and specular reflectivities. 

A geometric configuration is defined by the cus- 
tomer, consisting of a set of defining primary surfaces 
and their nodes. A primary surface may be a paral- 
lelogram, disc, cylinder, sphere, or cone. Nodes are 
defined as subsurfaces on the primary surfaces. One 
input data case consists of a geometric configura- 
tion together with commands for finding Script F’s 
between objects within the configuration. Commands 
are executed as they are encountered; when a Script 
F command is encountered, the following operations 
are performed: The node summation areas are cleared 
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to zero. These areas consist of “Script-F^ sums for 
1 and 2 wavelength nodes along with force and torque 
sums., The presence of these areas is entirely depend- 
ent upbn ^ of node being encountered. The 
counter for the total number of rays considered is set 
to zero. The counter is incremented by 1 . 

Notes: 

1. The program is written in FORTRAN IV language 
for use on the IBM 360 computer. 

2. At most, no more than 100 primary surface tests 
may be made per second. In addition, there is a 
relatively high machine-time cost for high-precis- 
ion parametric results. 

3. An advantage is the great generality of problems 
treatable and rapidity of solution from problem 
conception to receipt of results. 


4. 


Inquiries concerning this program may be directed 


Computer Center 
University of Georgia 
Athens, Georgia 30601 
Reference: B69- 10038 


Patent status : 

No patent action is contemplated by NASA. 

Source: James H. Scates and Jerry L. McWilliams 
of The Boeing Company 
under contract to 
Marshall Space Flight Center 
(MFS-15051) 
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